Disorder effects in the quantum kagome antiferromagnet ZnCu3(OH)6Cl2 



00 

o 
o 

(N 
>>■ 



(N 
(N 



M. J. Rozenberg^^^ and R. Chitra^ 

^Laboratoire de Physique des Solides, CNRS-UMR8502, Universite de Paris-Sud, Orsay 91405, France. 
^ Departamento de Fisica, FCEN, Universidad de Buenos Aires, 
Ciudad Universitaria Pah. I, (1428) Buenos Aires, Argentina. 
^ Lahoratoire de Physique Theorique de la Matiere Condensee, UMR 7600, 
Universite de Pierre et Marie Curie, Jussieu, Paris-75005, France. 

(Dated:) 

Motivated by the recent NMR experiments on ZnCu3(OH)6Cl2, we study the effect of non- 
magnetic defects on the antiferromagnetic spin-i kagome lattice. We use exact diagonahzation 
methods to study the effect of two such defects on finite size systems. Our results, obtained without 
adjustable parameters, are in good quantitative agreement with recent ^''O NMR data. They pro- 
vide support for the experimental interpretation of the presence of defects within the kagome layers 
due to Zn/Cu substitutions. Our results also show that disorder effects become relevant at lower 
temperatures, raising questions about the experimental evidence for the absence of an intrinsic spin 
gap in the kagome 2D layers. 

PACS numbers: 75.10.Jm, 75.40.Gb 
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The quest for spin liquids in frustrated quantum mag- 
nets remains one of the main challenges in condensed 
matter physics. Spin liquids are quantum spin systems 
where quantum fluctuations thwart any kind of ordering 
of the spins even in the limit of T — > 0. Typical criteria 
for a candidate system include: a high degree of frustra- 
tion, low dimensionality and small spin. All three fea- 
tures are realized in the spin-i Heisenberg antiferromag- 
net on a kagome lattice, which is a 2D network of corner 
sharing triangles. In this context, the recent discovery of 
the compound herbertsmithite, ZnCu3(OH)6Cl2 0, [3], 
presumed to be a perfect physical realization of such a 
model, has triggered a spate of experimental studies de- 
voted to unveiling the intrinsic behavior of the kagome 
antiferromagnet. A host of experimental techniques, in- 
cluding neutron scattering [sl, NMR[i, @| and ^SR^ 
have already been used to study the magnetic behav- 
ior of this compound. However, comparisons of experi- 
mental and theoretical results are hampered by the fact 
that not much is known about the kagome system ex- 
cept that the ground state is a singlet, the spectrum is 
characterized by a macroscopic quasi degeneracy of low 
lying singlet states, a very short correlation length and 
possibly a small gap to a macroscopic number of quasi 
degenerate triplet excitations 0, [1] ■ 

Surprisingly, the plethora of experimental data ob- 
tained show that magnetic degrees of freedom remain 
active at low frequencies and temperatures and these do 
not conform to the naive picture of a gapped system. 
This has generated a debate as to whether the experimen- 
tal systems are pure enough to reveal intrinsic kagome 
behavior, or, alternatively, whether additional interac- 
tions that go beyond simple nearest neighbor antiferro- 
magnetic coupling should be included. From a chemical 
perspective, the most delicate aspect is the control over 
the positions that the Cu and Zn atoms occupy. The 



former are magnetic and should be ideally confined to 
the 2D kagome planes, while the latter, non-magnetic, 
should occupy interlayer sites that effectively decouple 
the planes. For each exchange of these atoms, two defects 
are created: (i) the Cu that occupies the Zn site at the in- 
terlayer, remains rather weakly coupled and contributes 
an effective "free" spin-^ term to the bulk susceptibility, 
and (ii) the Zn that occupies a Cu site in a layer, creates a 
non-magnetic defect within the kagome plane which will 
likely modify the intrinsic kagome behavior. 

On the theoretical front, different ideas were proposed 
to account for the unexpected observations: in Ref. 
it was shown that the inclusion of Dyazolinshki-Moriya 
(DM) interactions may enhance the uniform susceptibil- 
ity at intermediate temperatures. On the other hand, 
other studies have shown that magnetic defects may also 
account for several experimental observations including 
the uniform susceptibility, the specific heat [lo| and the 
dynamical susceptibility 11]. 

To elucidate this situation, a rather powerful experi- 
mental technique is the NMR using oxygen isotopes. The 
oxygen atoms only occupy sites within the 2D planes and 
since the hyperfine coupling is very short ranged, oxygen 
NMR probes the physics of the plane and is insensitive 
to the spin-i magnetic contributions from misplaced Cu 
that occupy interlayer sites. Moreover, NMR probes both 
static and dynamic behavior. This decisive experiment 
was recently performed by Olariu et al. [l^ . The main 
results reported were (i) the existence of non-magnetic 
defects within the 2D layers, and (ii) the apparent ab- 
sence of a spin gap. To comprehend the effect of defects, 
here we study a spin-i antiferromagnetic kagome lattice 
model with non-magnetic impurities and obtain the theo- 
retical predictions for the temperature dependence of the 
NMR spectra and the nuclear spin relaxation rate 1/Ti. 
The quantitative agreement between our results and the 
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experimental data unambiguously demonstrates the ex- 
istence and importance of non-magnetic defects within 
the kagome layers of the herbertsmithite. However, we 
also find that the disorder due to these defects produces a 
smearing of the NMR line that is most significant at low 
temperatures, raising questions about the experimental 
evidence for the absence of a spin gap in the intrinsic 
kagome lattice system. 

The pure hamiltonian is given by the Heisenberg model 
on the kagome lattice. 



TC = J ^ ^ Si ■ Sj 



(1) 



where, S represent the spin-i operators and J the anti- 
ferromagnetic interactions between two nearest neighbor 
spins. Non-magnetic defects are introduced by eliminat- 
ing a spin operator at a given site. This model is not very 
amenable to analytic studies and is usually studied using 
numerical methods Here, due to the defects, we can- 
not take advantage of the symmetries of the pure lattice 
and moreover, as we are interested in the full tempera- 
ture dependence, we employ exact diagonalization. We 
therefore, consider a finite size system of 12 sites as the 
one schematically depicted in Fig. [1] with periodic bound- 
ary conditions. Note that the 12 site cluster (in the pure 
case) preserves all the original symmetries of the full lat- 
tice. In addition, since the magnetic correlation length 
in the pure system is extremely short, of the order of 
one lattice spacing [13], we expect that the 12 site clus- 
ter calculation may be adequate for the specific issues 
addressed in the present study 

We now study the effect of defects in the 12 site cluster. 
We do not consider the case of a single defect, since this 
would induce local moments which would then contribute 
a Curie-like term to the shift. No such contribution has 
been seen in the NMR study of Ref.[l3 which probes 
only the magnetic response in the kagome layers. In con- 
trast, a free spin-i contribution was in fact seen in a 
previous NMR study [1], using Cu and CI, which probed 
both in and out of the layer magnetic sites. Thus, the 
Curie-like contribution should be attributed to an inter- 
layer Zn/Cu substitution. In what follows, we consider 
all possible configurations of two non-magnetic defects. 
Fig. [T] shows an example of one such configuration. Two 
impurities in the twelve site cluster correspond to a defect 
concentration of 16% which is larger than the experimen- 
tal estimate of 6-10 % . However, due to the very short 
correlation length in the pure Kagome system, we expect 
our results to remain valid for a relatively large range 
of impurity concentrations, including the experimental 
estimate. 

The specificity of the ^^O NMR stems from the struc- 
ture of ZnCu3(OH)6Cl2, where oxygen is confined to the 
kagome planes. In a perfectly pure structure, each oxy- 
gen atom would have two equivalent magnetic Cu(2-|-) 
neighbors that are also located within the 2D kagome 




FIG. 1: (Color online) Schematic representation of the 12 
site kagome lattice with 2 non-magnetic defects. The bond A 
denotes a bond with a single spin and bond B denotes a bond 
with two spins. For a given configuration of defects, there are 
many inequivalent A-bonds and B-bonds. 



layer, and with whom it interacts via the hyperfine cou- 
pling. In an NMR experiment, a strong magnetic field is 
applied producing a local magnetization of the Cu ions. 
This magnetization modifies the strength of the hyper- 
fine coupling with the nuclear moment, producing a shift 
in the nuclear magnetic resonance line. Thus, the shift 
in the oxygen line is a direct measure of the individual 
magnetizations, and hence of the spin susceptibility, of 
the neigboring Cu sites (more precisely of the combined 
effect of the two neighbors). A knowledge of the local 
fields then gives us information about the existence of any 
underlying spin order. For example, if two neighboring 
spin sites were strongly antiferromagnetically correlated 
they would produce no net shift to the central oxygen 
line. In a pure system, since all oxygen sites are equiva- 
lent, the NMR spectra should have a single well defined 
peak with a temperature dependent shift. However, if 
non-magnetic defects are present, then the oxygen ions 
may encounter two types of magnetic environments: one 
comprising two nearest-neighbor Cu, and the other with 
only one Cu. (The case of two non-magnetic neighbors 
would produce no NMR shift.) This would then be di- 
rectly refiected in the NMR spectra as a superposition of 
two peaks, each one related to the respective contribu- 
tion arising from the two kinds of magnetic environments 
that an oxygen moment may see. 

We now present our theoretical predictions for the 
NMR. The only parameter in our model is the magni- 
tude of the antiferromagnetic coupling J, which we set 
to the experimentally estimated value ~ 170K. The lon- 
gitudinal magnetic field in our work H = 0.05J is of 
similar magnitude to the value of 7Tesla used in Ref. Il2. 
To obtain the shifts, we first compute the finite tem- 
perature magnetization at every site rrii. The magnetic 



shifts Kij at the "oxygen sites" (which in the physical 
structure are equidistant to two neighboring vertex of 
the kagome lattice) are then estimated as the sum of the 
two contributions rui and rrij of the induced moments at 
the two sites i and j of each bond of the lattice. The 
defect free 2-spin bonds, such as bond B, depicted in 
Fig. [T] have in fact two contributions; however, the 1- 
spin bonds, such as bond A, have a contribution from 



a sole magnetic moment rrii. Note that for any given 
configuration of defects, since the translational symme- 
try is broken, there are many inequivalent 1-spin bonds 
and 2-spin bonds. The 24 shifts Kij for each of the four 
non-equivalent two-defect configurations are computed, 
and the data at each temperature are condensed into a 
histogram for the distribution of the magnitudes of the 
KijS. These histograms can be directly compared to the 
NMR spectra obtained at different temperatures. [3] 

Fig. [5] summarizes our results for the predicted NMR 
spectra as a function of temperature. In a pure, defect 
free system, all sites would be equivalent and the his- 
tograms would simply show a single sharp peak, whose 
position depends on the temperature. Since the shift 
is proportional to the local magnetizations due to the 
external field, its magnitude is directly proportional to 
the uniform susceptibility. However, when defects are 
introduced, translational invariance is broken and a dis- 
tribution of local magnetic moments are induced, which 
results in the the histograms of Fig. [51 The main features 
that emerge from the data is the presence of two promi- 
nent "ridges" that can be inferred from the temperature 
dependence of the histograms. The is a clear manifes- 
tation of the presence of two types of magnetic environ- 
ments: some bonds have two spins (ie, no defect), while 
others have only one spin (and one defect). The position, 
or shift, of the rightmost ridge is larger by a factor of two 
compared to the other one. In fact, the former originates 
from contributions of bonds with two spins while the lat- 
ter from that of bonds with one spin and one defect. At 
high temperatures, as the spins are effectively decoupled, 
this scaling feature is intuitively expected. However, the 
fact that this scaling remains valid well below T/J ^ 1, 
is a non-trivial signature of strong frustration and an ex- 
tremely short correlation length, that are to be expected 
in a spin fluid state. Significantly, this behavior is in 
very good agreement with the experimental -'^^O NMR 
data |12| . 

In order to visualize the results under a different light, 
we redisplay the temperature dependent histograms as 
an intensity plot in Fig. [H For reference, we also show 
in the inset the calculated data for the uniform spin sus- 
ceptibility xC^) for the pure system with no defects. At 
higher temperatures, x shows the Curie tail, while at low 
T it shows an exponentially activated behavior, due to 
the presence of a gap in the energy spectrum. The main 
panel shows that the maximum intensity lincshapes of the 
two ridges roughly follow the behavior of x(T). However, 
in contrast to the x{T) of the pure system, the T depen- 
dence of both lineshapes become flat between 70 to IQQK 
where they attain their respective maximum values. Sig- 
nificantly, a similar flattening in the same temperature 
regime around T « 0.5J was experimentally observed in 
RefE 

Another interesting feature that our results of Figs. 
and [3] show, is the dramatic smearing of the peaks in 
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FIG. 2: Histograms for the values of the sum of the local 
magnetization of the spins on each lattice bond. The mag- 
netizations are proportional to the shifts of NMR. The his- 
tograms correspond to the theoretical estimate for the NMR 
spectra. 



the histogram at low temperatures. Nevertheless, be- 
low T w 5QK{Ri 0.3 J) one observes that while the lower 
ridge remains relatively well defined, the upper one shows 
a substantial dispersion. This behavior can be qualita- 
tively interpreted as follows: the lower ridge arises from 
A bonds which all share the common feature that their 
spins are adjacent to a defect site. In contrast, the higher 
ridge receives contributions from B bonds whose spins 
may be located at different distances from the defects 
and hence results in a significant smearing of the his- 
togram. This feature is also qualitatively observed in 
the experimental NMR data of Ref. 1 121. T here, the peak 
of the lower shift (denoted D in Ref. | l2| ) remains very 
sharp, while that of the higher shift (denoted M) becomes 
significantly rounded beneath 85K. This defect induced 
smearing, which is most significant a lower temperatures, 
casts doubts on the experimental evidence for the closing 
of the intrinsic spin gap in the kagome system 1^. In 
fact, in the measured NMR spectra, the peaks can no 



longer be resolved beneath lOK. Our results rather sug- 
gest a scenario where the apparent low lying magnetic ex- 
citations reported in the experiment should be attributed 
to disorder effects stemming from uncontrolled substitu- 
tion of Zn and Cu atoms in the structure of the kagome 
planes. 17] 

Finally, we study the NMR spin relaxation rate esti- 
mated as 



^-1 ^ lmx(c^o) 



(2) 



where ujo is the Larmor frequency of the applied external 
magnetic field and x{^) is the dynamical susceptibility. 



FIG. 3: (Color Online) Intensity plot of the histograms of 
Fig. (2] K denotes the theoretical estimates for the shifts (in 
arbitrary units). The inset shows the uniform magnetic sus- 
ceptibility x{T) for the pure 12 site system. 



We compute x(w) for all sites and all defect configura- 
tions and evaluate the average value of 1/Ti. We should 
note, however, that unlike the magnetization previously 
computed for the NMR shifts, the NMR rate probes the 
low frequency behavior of the dynamic susceptibility and 
consequently more prone to finite size effects. Nonethe- 
less, we expect our results to provide useful qualitative 
insight. Our calculations for the finite temperature relax- 
ation rate are shown in Fig. [H We plot the data for four 
particular defect configurations denoted by the Manhat- 
tan distance between the two defects ( Ml, M2, M3 and 
M3' ) and also the average of the four sets of data (mul- 
tiplied by 4 for easier visualization). In the inset of the 
figure, we show the results as a function of inverse tem- 
perature 1/T on a semi-log scale. The data for Ml, M2, 
M3 and MS' are linear at low T , which is expected from 
simple activated behavior. In contrast, the average data 
(solid red line) clearly deviates from simple activated be- 
havior, displaying a gentle upward curvature. These re- 
sults are consistent with the experimental data of Ref.ll2. 
The lack of linearity of the measured 1 /Ti and a putative 
power law fit to the data was interpreted as a signature of 
the closing of the intrinsic spin singlet/triplet gap in the 
system. However, our results indicate that an alterna- 
tive scenario is that the non-linear behavior is associated 
with the existence of a multitude of small activation gaps 
resulting from the presence of defects. 

To conclude, we have studied the effects of non- 
magnetic impurity defects on a finite size kagome lat- 
tice using exact diagonalization techniques. With no ad- 
justable parameters our results show a remarkable agree- 
ment with recent ^^O NMR data and validate the exper- 
imental interpretation regarding the origin of two promi- 
nent lineshifts in the NMR spectra that were associated 
with Zn/Cu substitution in the 2D kagome planes. In 
particular, the model results reproduce the factor of two 
scaling between the lineshifts which persists down to low 



FIG. 4: (Color online) The relaxation rate l/Ti as a func- 
tion of T. Ml, M2, M3 and M3' denote configurations with 
the defects at distances 1, 2, 3 and 3 respectively (there are 
two inequivalent realizations for distance 3). The solid line is 
the sum of the four contributions. The inset shows the same 
quantities on a semi-log scale as a function of inverse temper- 
ature. Straight lines indicate simple activated behavior. The 
dashed (blue) straight line is drawn for reference. 



temperatures which is indicative of the spin liquid nature 
of the underlying system. It also captures the peculiar 
temperature dependence of the measured spin suscepti- 
bility and provide insight on the unequal smearing of the 
two lineshifts. On the other hand, our results also cast 
important doubts on the experimental evidence for the 
absence of a spin gap in the system, suggesting that dis- 
order strongly affects the low temperature behavior. This 
calls for renewed efforts to achieve a better chemical con- 
trol on sample quality. 
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